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Abstract: We present an implementation for slepton pair production at hadron colliders 



in the POWHEG BOX, a framework for combining next-to-leading order QCD calculations 



with parton-shower Monte-Carlo programs. Our code provides a SUSY Les Houches Ac- 
cord interface for setting the supersymmetric input parameters. Decays of the sleptons 
and parton-shower effects are simulated with PYTHIA. Focussing on a representative point 
in the supersymmetric parameter space we show results for kinematic distributions that 
can be observed experimentally. While next-to-leading order QCD corrections are sizable 
for all distributions, the parton shower affects the color- neutral particles only marginally. 
Pronounced parton-shower effects are found for jet distributions. 



1 Introduction 



With the start-up of the CERN Large Hadron Collider (LHC) unprecedented opportunities 
have emerged for experimentally accessing the terascale. The capability of the LHC was 
impressively demonstrated by the recently announced discovery [1, 2] of a new particle 
compatible with the Higgs boson of the Standard Model. Yet, numerous extensions of the 
Standard Model are being discussed, some of the best-motivated ones being supersymmetric 
(SUSY) models that predict the existence of supersymmetric partners for all Standard 
Model particles. 

Direct production of the scalar partners of the leptons proceeds via slepton pair pro- 
duction, if R-parity is conserved. Compared to the production of the color-charged squarks 
and gluinos, the production cross sections are smaller but also the signatures in the de- 
tector are cleaner (see, e.g., Ref. [3] for a recent review). Current searches at ATLAS [4] 
already supplement the LEP limits [5] by excluding masses of the left-handed sleptons up 
to 185 GeV for exclusive decays into a lepton and the lightest neutralino, depending on 
the neutralino mass. Precision measurements in the slepton sector would be possible at a 
future e + e~ linear collider [6, 7] and to a less general extent also at a muon collider [9]. 

For slepton pair production at hadron colliders, the leading order (LO) [8, 10] and 
next-to-leading order (NLO) QCD [11] and SUSY-QCD [12] contributions are known for 
the total cross sections as well as for differential distributions. The latter have been imple- 
mented in the computer package PRDSPIND [13], which is publicly available. Uncertainties 
inherent to the fixed-order parton-level calculations have been further reduced by taking 
resummation effects into account. In Refs. [14] and [15], transverse momentum and thresh- 
old resummation for slepton pair production at hadron colliders was considered, while 
Ref. [16] provided joint resummation predictions for this class of reactions. More recently, 
using methods of soft-collinear effective theory invariant-mass distributions and total cross 
sections have been presented at next-to-next-to-next-to-leading logarithmic accuracy [17]. 

Ideally, precise predictions for the production cross sections should be combined with 
realistic simulations of parton-shower effects, multi-parton interactions, and underlying 
event - effects that are omni-present at hadron colliders. Moreover, the complex decay pat- 
terns of the heavy SUSY particles should fully be taken into account and realistic analysis 
cuts should be applied. These requirements call for the use of multi-purpose Monte-Carlo 
generators such as HERWIG [18, 19] or PYTHIA [20], which are, however, mostly restricted to 
leading order matrix elements for hard scattering processes. A framework that allows to 
combine the merits of both, flexible parton-shower programs with precise NLO-QCD cal- 
culations, is the so-called POWHEG approach [21, 22]. The POWHEG BOX [23] is a public tool 
providing all general building blocks of this method, but requiring the user to individually 
implement process-specific pieces such as matrix elements for the hard scattering process 
at NLO-QCD accuracy. Several Standard-Model processes bearing similar kinematic fea- 
tures as slepton pair production are available in the POWHEG BOX [24, 25]. Selected SUSY 
processes have been studied in the POWHEG approach as well [26, 27]. 

In this article, we present an implementation of charged slepton pair production at 
hadron colliders in the POWHEG BOX and study the relevance of parton-shower effects for 
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representative distributions. The slepton decay implementation of PYTHIA is used to pro- 
duce experimentally accessible final states. Firstly, we describe our calculation in Sec. 2. 
In Sec. 3, we provide representative phenomenological results for a specific parameter point 
of the minimal supersymmetric extension of the Standard Model (MSSM). Our conclusions 
are given in Sec. 4. 

2 Framework of the calculation 

The NLO-QCD and SUSY-QCD (SQCD) corrections to slepton pair production at hadron 
colliders are known for quite some time [11, 12] and publicly available, e. g. in the parton- 
level Monte Carlo program PROSPINO [13]. Nonetheless, we decided to re-calculate the 
full 0(a s ) corrections to this class of processes in the framework of the MSSM. In the 
following, we consider the production of a pair of charged sleptons and leave a discussion 
of final states involving sneutrinos for future work. 

At LO, the production of charged slepton pairs at hadron colliders dominantly proceeds 
via the Drell-Yan type annihilation of a quark-antiquark pair into a virtual Z boson or 
photon that in turn decays into a pair of heavy sfermions, c.f. Fig. 1 (a). The SQCD- 
NLO corrections to this reaction comprise virtual one-loop corrections and real-emission 
contributions with an extra parton in the final state. To the former only self-energy and 
vertex corrections contribute, since the sleptons do not carry color charge, see Fig. 1 (b). 
For the real-emission corrections we consider qq annihilation diagrams with a gluon being 
radiated off either of the incoming fermions, and crossing-related contributions where the 
gluon is promoted to the initial state while the final state features a quark or anti-quark 
in addition to the slepton pair, as shown in Fig. 1 (c). 

We employ the on-shell scheme for the renormalization of the wave functions of the 
massless quarks, which requires the evaluation of additional self energy diagrams not de- 
picted here. Other quantities such as the couplings receive no non-trivial contributions at 
the order of perturbation theory considered. We account for large universal corrections 
to the electromagnetic coupling due to light fermions by a running coupling. We input 
its value at the scale mz and use the resummed one-loop contributions to evolve to the 
invariant mass of the slepton pair. 

We generate the Feynman diagrams with QGRAF [28] and build the interference terms 
with in-house developed FORM [29, 30] scripts. The relative sign of interference produced 
by QGRAF cannot be used due to the presence of Majorana fermions; we generate this 
sign according to Ref. [31]. Contraction of Lorentz indices leaves us with scalar integrals 
which are reduced to master integrals with the program Reduze2 [32, 33]. The resulting 
expressions are exported to Fortran 77 files that make use of the QCDLoop library [34] 
for the numerical evaluation of the scalar master integrals. Conventional dimensional 
regularization is used to regularize both ultraviolet (UV) and infrared (IR) divergences. 
No supersymmetry-restoring counterterms are needed here [35]. The calculation employs 
7 5 = («/4!)e A1I/po -7 M 7 l/ 7 p 7 " for the axial couplings [36, 37]. Additionally we performed an 
independent calculation utilizing the Mathematica package FeynArts [38] for amplitude 
generation, FormCalc [39, 40] for algebraic simplifications including tensor reductions, and 
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Figure 1. Fcynman diagrams contributing to slepton pair production at hadron colliders: tree-level 
contribution (a), virtual NLO-QCD and SQCD corrections (b) and real-emission contributions (c) 
from different channels. 

LoopTools [39, 41] for the numerical evaluation of the scalar integrals. This calculation uses 
the so-called naive anticommuting scheme [42] for 7 5 . The results of these two calculations 
are in complete agreement with each other. In both approaches, we neglect masses of 
quarks and leptons for the first and second generations. Consequently, the scalar partners 
of these left and right-chiral fermions are treated as mass eigenstates. While partons of 
the third generation may safely be ignored in the initial state, mixing of staus is taken into 
account. 

Having calculated the 0{a s ) corrections to the partonic scattering process has put us 
in a position to work out an implementation of slepton pair production in the POWHEG BOX. 
We merely have to provide the following ingredients: 

• a list of all flavor structures for the LO and the real-emission contributions, 

• a suitable parameterization of the Born phase space, 

• the Born amplitudes squared for all partonic subprocesses, 

• the color correlated and the spin correlated Born amplitudes, 

• the finite part of the virtual QCD and SQCD corrections, 

• the real-emission matrix elements squared for all partonic subprocesses. 

Once these building blocks are implemented in the POWHEG BOX, the program itself performs 
the phase-space integration and convolution with the set of parton-distribution functions 
selected by the user. IR divergent configurations are taken care of internally by a subtrac- 
tion procedure based on Ref. [43]. The POWHEG BOX also provides the means for checking 
that the real-emission contributions approach the respective counterterms for soft and 
collinear configurations. This constitutes a useful debugging tool as it tests the correct 
normalization of the real-emission amplitudes relative to the Born matrix elements that 
are used for the computation of the counterterms. 
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We have compared our total cross section with results obtained from PROSPINO 2.1 
and found agreement both at LO and at NLO. In addition, we checked that our virtual 
Standard-Model-QCD and SUSY-QCD corrections as well as LO and NLO total cross 
sections are in agreement with [17]. 

In our implementation, spectrum dependent decays of the on-shell sleptons are incor- 
porated via PYTHIA 6.4.25. For the input of fully general MSSM parameters, we employ 
the SUSY Les Houches Accord (SLHA) [44, 45] interface. The production part of our im- 
plementation utilizes the SLHALib2 library [46] (SLHA 1 & SLHA 2), while the decay code 
of PYTHIA already features configuration via SLHA 1. In this way, the output of spectrum 
generators can be passed directly to our program in a well-defined and user-friendly way. 

3 Phenomenological results and discussion 

Our implementation of slepton pair production in the POWHEG BOX will be made publicly 
available at the homepage of the POWHEG BOX project, http://powhegbox.mib.infn.it. 
We encourage the phenomenologist to use this code with his or her own preferred settings for 
the SUSY spectrum, Standard-Model input parameters and experimental selection cuts. 
A documentation of the code including recommended values for technical parameters is 
provided together with the code. Here, we show representative results for a specific point 
in the SUSY parameter space. 

We consider proton-proton collisions at the LHC with a center-of-mass energy of \/~S = 
8 TeV. For the electroweak parameters we use input provided by the particle data group 
(PDG) [47]. We set the mass of the Z boson mz = 91.1876 GeV, the fine-structure 
constant a = 1/137.036 and use a running electromagnetic coupling in the on-shell scheme 
with a(mz) = 1/128.919. For the electroweak mixing angle we obtain sin 2 9w = 0.23103 
from the PDG value for Fermi's constant, mz and a{mz) using the tree-level relation. The 
Z width is set to its PDG value but has no impact on the results discussed in the following. 
For the masses of the sleptons and the lightest neutralino we set 

m l a = 180 GeV for Ir = eR,p,R, 

m -o = 80 GeV, (3.1) 

while for squarks and gluino we use m,q L = m g - fl = 1500 GeV for q = u,d,c,s and m g = 
2000 GeV. While SUSY particles of these masses are not excluded by collider data, from 
Ref. [4] one may expect such a slepton sector to be accessible by future slepton searches at 
the LHC experiments. The virtual SQCD corrections induced by the squarks and gluinos 
are included but found to be numerically irrelevant because of their large masses in this 
setup. 

For the numerical analysis we restrict ourselves to the pair production of the SUSY 
partners of the right-handed leptons of the first and second generation, the R-selectrons 
and R-smuons. The sleptons are assumed to exclusively decay into a lepton and the light- 
est neutralino, where we employ PYTHIA 6.4.25 for the decay. Throughout, we switch off 
QED radiation, underlying event and hadronization effects in PYTHIA. Unless stated other- 
wise, factorization and renormalization scales are identified with the invariant mass of the 
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Figure 2. Invariant- mass distribution of the slcpton pair (a) and transverse momentum of the 
negatively charged selectron (b) at LO (green), NLO-QCD (black), and with POWHEG+PYTHIA (red) 
for our default LHC setup with VS = 8 TeV. 



slepton pair. For the parton-distribution functions of the proton, we use the MSTW2008 
parameterization [48] as implemented in the LHADPDF library [49] . The real-emission contri- 
butions of the NLO-QCD calculation as well as the parton shower can give rise to partons 
in the final state. These are recombined into jets according to the anti-fc^ algorithm [50] 
as implemented in the FAST JET package [51, 52], with R = 0.4 and y iCt < 4.5. 

For the total cross section we find a = 5.93 fb at NLO-QCD. In order to illustrate 
the dependence of our results on unphysical scales, we evaluated the total cross sections 
for our default setup at LO and NLO-QCD, varying the factorization and renormalization 
scales in the range /j,q/2 to 2hq, with hq = ^ m i R - Note that the renormalization scale does 
not enter the LO results, since at tree-level slepton pair production is a purely electroweak 
process. While at LO we find a scale uncertainty of about 8%, the NLO results for different 
values of fio differ by only 4%. 

In Figs. 2 (a) and (b) we present the invariant-mass distribution of the slepton pair 
and the transverse momentum of the negatively charged selectron, respectively, for our 
default setup at LO, NLO, and for POWHEG+PYTHIA (NLO+PS). The statistical errors on 
the distributions are at the order of 1% per bin or smaller and therefore not explicitly 
indicated in the histograms. As expected, relatively large positive NLO-QCD corrections 
are encountered for both observables, while the impact of the parton shower on the NLO 
result is marginal for these slepton distributions. 

More pronounced parton-shower effects can be observed in distributions related to the 
hardest jet produced in association with the slepton pair. In the NLO-QCD calculation, 
a jet can only arise from the final-state parton of the real-emission contributions. In 
contrast, the parton shower can give rise to events with even more than one hard jet in 
the POWHEG+PYTHIA result. The transverse-momentum distribution of the hardest jet is 
shown in Fig. 3 (a) for the fixed-order and the POWHEG+PYTHIA calculation. While in the 
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Figure 3. Transverse momentum (a) and rapidity distribution (b) of the hardest jet for slepton 
pair production at the LHC with y/S = 8 TeV. In (b) we additionally require p^f* > 20 GeV. 




Figure 4. Transverse-momentum distribution of the negatively charged lepton (a) and rriT2 dis- 
tribution (b) with POWHEG+PYTHIA within the cuts of Eq. (3.2) for slepton pair production at the 
LHC with VS = 8 TeV. 

fixed-order QCD prediction eta/dp!^ rises steadily towards very low transverse momenta, 
this rise is damped by the Sudakov factor in the POWHEG+PYTHIA result. In either case, the 
extra jet tends to be produced at central rapidities, as illustrated by Fig. 3 (b). There, 
we imposed an extra cut on the transverse momentum of the jet, pj? > 20 GeV to avoid 
contributions from very soft partons produced by the parton shower. 

An important feature of our code is the possibility to interface the calculation of the 
slepton pair production cross section with PYTHIA in such a way that decays of the sleptons 
can be simulated, and the user gets access on kinematic distributions of the decay products 
within selection cuts of his own choice. In Fig. 4 (a) we present predictions at the NLO+PS 
level for the transverse-momentum distribution of the hardest negatively charged lepton 
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within our standard setup, and with additional cuts on the charged leptons, requiring them 
to be hard and centrally produced, 

j4 > 20 GeV and rf < 2.5 . (3.2) 

Providing predictions in terms of final decay products, our code can be used directly also 
for more elaborate analysis techniques such as studies of the distribution defined in 
Ref. [53]. This is demonstrated in Fig. 4 (b), where we used the code of Ref. [54] for the 
calculation of m^- 

While these results are meant as an illustration of the capability of our code, we would 
like to stress that the user is free to add distributions, modify cuts, and set masses of the 
SUSY particles according to his needs in a straightforward manner. 



4 Conclusions 



In this work, we have presented an implementation of slepton pair production at hadron col- 
liders in the POWHEG BOX, a framework that allows to combine NLO-QCD calculations with 
parton-shower Monte-Carlo programs. We have described how we calculated the QCD and 
SQCD corrections to the hard scattering process, and then explained the technical details 
of their implementation in a publicly available computer code. The program we developed 
contains an SLHA interface that allows the use of MSSM input parameters provided by 
stand-alone spectrum generators. When our program is run with PYTHIA, decays of the 
sleptons can be simulated by the Monte-Carlo generator. 

We provided representative results for various kinematic distributions of the sleptons 
and their decay products for a specific point in the MSSM parameter space. Our analysis 
showed that NLO-QCD corrections modify tree-level predictions for slepton pair production 
processes significantly. Parton-shower effects can be pronounced for jet observables, while 
they are typically small for distributions related to the sleptons. 



Note added in proof 

During the completion of this work, a similar study appeared [55] , containing a simulation 
of slepton pair production at NLO-QCD in HERWIG++ [56]. Two separate implementations 
of this reaction in completely different program packages provide excellent means to check 
the independence of predictions from a specific tool. In this sense, it would be interesting 
to compare the results of the two implementations in the future. 
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